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We study gravitational collapse of the axion/dilaton field in classical low en- 
ergy string theory at the threshold for black hole formation. A new critical 
solution is derived that is spherically symmetric and continuously self-similar. 
The universal scaling and echoing behavior discovered by Choptuik in gravi- 
tational collapse appear in a somewhat different form. In particular, echoing 
takes the form of 5L(2,R) rotations (c/. ^-duality) . The collapse leaves be- 
hind an outgoing pulse of axion/dilaton radiation, with nearly but not exactly 
flat spacetime within it. 
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The striking numerical results of Choptuik |l] on the spherically symmetric gravitational 
collapse of a real scalar field, 0, have inspired an upsurge of interest in gravitational collapse 
just at the threshold for formation of black holes. Further numerical results for vacuum rela- 
tivity with axial symmetry and perfect fluids |3|-[5|] suggest that the phenomena discovered 
by Choptuik are not just restricted to spherically symmetric real scalar fields. 

The thought experiment employed by Choptuik, in the context of numerical relativity, 
is to "tune" across the critical threshold in the space of initial conditions that separates 
the non-black hole endstate from the black hole endstate, and to carefully study the critical 
behavior of various quantities at this threshold |TJ]. Two distinct kinds of scaling appear. 
First, just above the black hole threshold, the black hole mass is M hh (p) oc (p— p*) 7 , 7 ~ 0.37, 
where p is some tuning parameter of the initial conditions. 

Second, exactly at the threshold, there appears a unique critical solution acting as an 
attractor for all nearby initial conditions on threshold. This critical solution — which we 
will call a "choptuon" — exhibits a striking, recurrent "echoing" behavior: Asymptotically 
at small time-scales t and length-scales r near the collapse, it repeats itself at ever-decreasing 
scales t' = e~ nA t, r' = e~ nA r (for n e Z + ). Here A is either a fixed constant of the solution 
(A m log 30 for a real scalar field A « log 5 for vacuum gravity ||) demonstrating a dis- 
crete self- similarity; or an arbitrary constant, demonstrating a continuous self-similarity 
The choptuon is itself a smooth solution of the field equations to the past of the asymptotic 
point (t, r) = (0,0), but clearly has some kind of spacetime singularity visible from infinity 
at that point; for instance, the Riemann curvature diverges at that point. 

Choptuons represent in principle a means by which effects of Planck scale gravity, 
even quantum gravity or superstring theory, might be observable in the present universe. 
Therefore, the quantum gravitational and stringy generalizations of the choptuon demand 
study Ufl. 

In this Letter we report on threshold behavior of gravitational collapse, and Choptuik 
scaling, in classical 3+1- dimensional low-energy effective string theory: i.e., in general rel- 
ativity coupled to a dilaton $ and an axion a. The fields a and $ can be combined into a 
single complex field r = a + ^e -2 *, and then the effective action for the model is (omitting 
gauge fieldsQ) 

S= f^^jR-l^) . (1) 



2 (Imr 

The equations of motion from this action are 

R ab ~ x 9 (d a rd b f + d a fd b r) = , (2a) 

W„ T+ ^ = 0. (2b) 
Imr 

Because the model (|2|) arises from an underlying theory of quantum gravity, it should be an 
excellent starting point for investigating quantum effects. 



One might think that the choptuon could carry gauge charge as well as mass, however this 
appears not to be so. 
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We assume spherical symmetry. We also assume continuous self-similarity, which means 
we work exactly at the threshold for black hole formation, and construct a critical solution, 
as in HH . Restricting to a = in fl2|) reduces to the model studied by Choptuik, which 
has discrete self-similarity. Generic initial conditions for a will presumably lead at threshold 
either to discrete or to continuous self-similarity; further numerical work is needed to decide. 

The model (§) has an extra global symmetry not present in general relativity, an SL(2, M) 
symmetry that acts on <3> and a, but leaves the spacetime metric invariant; this is a classical 
version of the conjectured SL(2, Z) symmetry of heterotic string theory called S'-duality ||. 
It acts on t as 

ar + b 

where (a, b, c, d) £ R with ad — be = 1, while leaving invariant. 

The main results of this Letter are: (1) Choptuik echoing can and does occur in this 
model, taking the novel form of these SX(2,K) rotations associated with S'-duality; (2) A 
critical solution exists that is continuously self-similar, as in |3|||; (3) The global structure 
of the critical spacetime is determined. 

By the assumption of continuous self-similarity, there exists a homothetic Killing vector 
£ that generates global scale transformations 

C^gab = 2g ab . (4) 

How should t transform under £? The naive answer (from dimensional analysis) is that r 
should be invariant. However, quite generally, one must allow for spacetime symmetries to 
mix with internal symmetries. For instance, the Maxwell equations for the potential are 
conformally invariant only in the sense that a gauge transformation be allowed to accompany 
each conformal transformation. In our case, we must allow for r to be scale invariant up to 
an infinitesimal SL(2, M) transformation 

£ 5 r = a + d\T + a 2 r 2 , (5) 

with a, G R, and this will be crucial as it allows echoing to occur. 
In spherical symmetry the metric can be taken as 

ds 2 = (1 + u(t, r)) (-b(t, rfdt 2 + dr 2 ) + rW . (6) 

The time coordinate is chosen so that gravitational collapse on the axis of spherical symmetry 
first occurs at t — 0, and the metric is regular for t < 0. This metric remains invariant in 
form under transformations t = t(t') of the time coordinate, and this invariance can be used 
to set b(t, 0) = 1 for t < on the axis. By regularity (no cone singularity), u(t, 0) = on the 
axis for t < as well. Initial conditions, assumed smooth, may be given on the hypersurface 
t = -1. 

In these coordinates £ = td t + rd r . Continuous self-similarity (P is then implemented by 
defining a scale invariant coordinate z = —r/t. In terms of z, 

b(t,r)—b(z), and u(t,r)—u(z), (7) 
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where the real functions b(z), u(z) are functions solely of the scale invariant coordinate z. 
From Eq. (|J), the r field must satisfy 

r i, r = i- . \) , 8 

where uo is a real constant, so far arbitrary, and where the complex function f(z) depends 
solely on z.f^ 

We can now use Eqs. in the equations of motion (0). After some work, u(z) drops 
out, giving the system of ODEs 

z(fr 2 -z 2 ) -, zo;(6 2 -* 2 ) -, - C^|/| 2 

6(1-|/| 2 ) 2/J 6(1- |/| 2 ) 2lJ/ //J 6(1 -|/| 2 ) 2 ' 1 J 

n _ r ^(6 2 + ^ 2 ) a -, 2 / ^(6 2 + ^ 2 ) \ ^(6 2 + 2^ 2 ) 

7 6 2 (1- |/| 2 ) 2/ 1 (1- |/| 2 ) V 26 2 (1 - |/| 2 ) y 11 6 2 (1 - |/| 2 ) 2// 7 

2 / ^ 2 (1 + |/| 2 ) ^V|/| 2 \ o; 2 z 2 -, 

T~ . 1 ' ..9\/l I XI9\ L9/7..9 ..9\/l I XI9N9 I J L9 /1 I I-\9\ r >J J 



(b 2 -z 2 )(l - |/| 2 ) 6 2 (6 2 -^ 2 )(l- I/I 2 ) 2 / J 6 2 (1-|/| 2 ) 2 

i ^ A Mi + l/l 2 ) ^ 2 I/I 2 f9bl 

+ (6 2 ~M (l-l/l 2 ) ^(l-l/l 2 ) 2 ;' 7 - 1 J 

This system apparently has five singular points, as follows. First, the points z = ±0 just 
represent the usual axis r = of spherical symmetry, and regularity is easily imposed. 

The point z = oo corresponds to the hypersurface t = 0. However, spacetime should 
be smooth on this hypersurface, except at the axis, since it lies in the Cauchy development 
of the initial conditions. Indeed, the system (^) can be rendered regular at z = oo by the 
following change of variables: 

dw = b(z) — , w = at z = oo , (10a) 

F(w) = z-^f(z) , (10b) 

v( w ) = , (10c) 

u{w)=u{z). (lOd) 

In terms of the new independent variable w, the equations of motion become 

° = V ' + (TWf F ' P ' + 1 " (l-'lFl 2 ) 2 ' (Ua) 

2vF> + iuF F , pi+ 2FF" 2iuv ( 1 + |F| 2 u _^_r_\ — | pi 

(1-|F| 2 ) 2 1-|F| 2 (t; 2 -l)(l ' ' 1 
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2 We have adopted a global <SX(2,R) gauge here. In principle, u> can also be purely imaginary or 
zero, but these do not seem to give interesting solutions. 



4 



Finally, there are singular points wherever b(z) = ±z (or v(z) = ±1), at which the 
homothetic killing vector becomes null; such characteristic hypersurfaces always occur in 
similarity solutions to hyperbolic equations, and in this case describe light cones in spacetime, 
at which new data could be injected in the r field. The singular point b(z + ) = z + (or 
v(w + ) = 1) describes the past light cone of the spacetime singularity at (t,r) = (0,0); see 
Fig. 1. Since this cone is within the Cauchy development of the spacetime initial conditions, 
the solution must be smooth across it. 




FIG. 1. Interpretation of the axion/dilaton choptuon. There are horizons at v = +1 (character- 
istic hypersurface) and v = —1 (Cauchy horizon). There is a spacetime singularity at (i = 0, r = 0) 
which is a single point in the causal structure. Region I is a collapsing sphere of gravitationally 
bound field; in Regions Ha and lib this blends smoothly into an outgoing wave. The outgoing 
wave oscillates infinitely many times approaching the Cauchy horizon at v = —1. The dashed lines 
indicate the first two oscillations. Region III is nearly but not precisely flat. 

The singular point b(zS) = —z_ or v(w_) = —1 describes the future light cone of the 
spacetime singularity; see Fig. 1. This cone is not within the Cauchy development of the 
spacetime initial conditions; in particular, it is the Cauchy horizon. Therefore the solution 
may not be smooth there, and we should not impose any boundary conditions beyond 
continuity of r. Observers on this cone will see data coming from the singularity, and indeed 
it is a very interesting question to ask what they will see. Evolution across a Cauchy horizon 
can never be unique, but we will show below that there exists a unique extension that is 
spherically symmetric, self-similar, and smooth on the future axis t > 0. We expect that a 
generic subcritical solution near the choptuon will be similar to the extension through z_, 
and a generic supercritical solution will form a singularity near z_ . 

The main advantage of the self-similarity hypothesis is that it reduces the equations of 
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motion @ to ordinary differential equations, @ and ([11]). Unfortunately, these ODE's are 
still too complicated to solve explicitly. Thus we use numerical methods to determine a 
solution. 

We are looking for a solution that is well defined everywhere, and smooth everywhere 
except perhaps at the Cauchy horizon z_. Define two real functions f m (z) and f a (z) such that 
f(z) = fm(z)e l f a ( z \ At any regular point z r of @, we need to provide the initial conditions 
b(z r ), f m (z r ), f^Zr), and fa(z r ) to begin numerical integration (Eq. (g) is independent of 
f a (z r )). At any singular point, we must take care to ensure that the fields (and their first 
derivatives) are continuous. As mentioned before, rescaling t allows us to set b(z = 0) = 1 
when t < 0. Requiring regularity at z = then fixes /^(0) and f' a (0) in terms of / m (0) = 
|/(0)|. A Taylor expansion of the fields b(z) and f(z) around z = then provides initial 
conditions for integrating outwards from z = €q, a small number. 

We can numerically integrate out from z = eo until the system reaches the next singular 
point at z + where b(z + ) = z + . The system must be smooth at z + , but this is numerically 
difficult to impose by just integrating out from z — e . So instead, we require that b and / be 
continuous at z + , and perform another Taylor expansion around z + , demanding smoothness. 
This provides initial conditions for integrating inwards from z = z + — e + in terms of |/(z + )| 
and z + . We then try to match the integration out from z = €q to the integration in from 
z = z + — e + at an intermediate point, which we choose to be z = 1. Each integration gives 
four final conditions, and there are four parameters to adjust, |/(0)|, u, z + , and so 
we have a well-determined system of equations. We find that if 

uj = 1.17695272200 ± 0.00000000270 , (12a) 

z + = 2.60909347510 ± 0.00000000216 , (12b) 

|/(0)| = 0.892555411872 ± 0.000000000224 , (12c) 

\f(z+)\ = 0.364210875022 ± 0.000000000760 , (12d) 

then the equations are smooth at z = and z = z + . 

The uncertainties quoted in fli~2]) represent one a numerical uncertainties. A good under- 
standing of the accuracy of the parameters is crucial especially near z_ where the spacetime 
turns out to be nearly, but not quite flat. To determine the errors in (|i~2"D , we have performed 
a large number of matchings using a powerful numerical integrator, starting each time with 
different initial estimates for the parameters, and with different values for e and the error 
tolerances for the integrator. All Taylor expansions are carried out to fourth order, and 
should not introduce any further uncertainties. 

The above equations uniquely determine the system at z = z + . We can therefore change 
to the w,v,F coordinate system (|T0|) and integrate from w = w + until the next singularity 
at w — w_ where v(w_) = — 1. This takes us through z = oo, t = into the region where 
t > 0. Near the Cauchy horizon at w_, F — > const, but F' oscillates rapidly. To see whether 
it is possible to continue through w_ , or whether it is a true singularity, we must look at the 
behavior of F in more detail. Near w_, eq. ([TTJb) becomes 



(w - w_)F" - c F' - ci = , (13) 

which has the solution 
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F(w) = c 2 + c 3 (w - w_) 1+c ° - —(w - w_) (14) 

Co 

where the are complex constants. Thus, F(w) will be continuous at w_ if Rec > — 1, 
and F'(w) will be continuous at w_ if Reco > 0. From eq. (PI), 

Rec = ^ 2 |i>-)| 2 = gfrjj ri5l 

(l-|F(w_)| 2 ) 2 -u; 2 |F(w_)| 2 1 -«(«;_)" 1 J 

Since \F(w_)\ is small, Rec > 0, and both F(w) and F'(w) are continuous at w_. Numer- 
ically, we find that 

w_-w + = 2.29540547918 ± 0.00000000192 , (16a) 
\F(w_)\ = 0.007603132483 ± 0.000000000125 , (16b) 
u (w_) = 0.0000800854142 ± 0.0000000000212 . (16c) 

Since u gives the mass aspect m = 2 a+u) ' ^ ne Cauchy horizon nearly, but not quite, carries 
data for flat spacetime. Because Reco is only slightly larger than zero, the function F(w) 
is just barely C 1 at w_. F" and all higher derivatives of F are discontinuous at w_. This is 
physically acceptable as the spacetime only depends on F and F' . 

The remaining region to consider is > z > z_. Since only / and /' are continuous at 
z = we can not repeat the same procedure that we used for < z < z + . Instead, we 
must begin integration at z = — eo and integrate out to z_. Since numerically we only know 
\f(z_)\, we can try to achieve that final condition in the integration. As before, rescaling t 
allows us to set 6(0_) = 1 and requiring regularity at z = 0_ fixes /^(0) and /„(0) in terms 
of |/(0_)|. This is a standard shooting problem with one free parameter and one target. We 
find that 

z_ = -1.000037262799 ± 0.000000000113 , (17a) 
|/(0_)| = 0.011742998497 ±0.000000000195 . (17b) 

This completes the last part of the spacetime. 

In conclusion, full scale numerical work on this model is highly desirable, for several 
reasons. The solution we have given is, by assumption, continuously self-similar. It is 
important to see if discretely self-similar solutions also exist — and if so which is the stronger 
attractor. Discretely self-similar solutions may manifest themselves by a form of Choptuik Q 
"echoing" in which the system is invariant after a finite scale transformation by some scale 
factor exp(— A), up to an SL(2,M) transformation. 
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